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How can one detect entanglement between multiple optical paths sharing a single photon? We 
address this question by proposing a scalable protocol, which only uses local measurements where 
single photon detection is combined with small displacement operations. The resulting entanglement 
witness does not require post-selection, nor assumptions about the photon number in each path. 
Furthermore, it guarantees that entanglement lies in a subspace with at most one photon per 
optical path and reveals genuinely multipartite entanglement. We demonstrate its scalability and 
resistance to loss by performing various experiments with two and three optical paths. We anticipate 
applications of our results for quantum network certihcation. 

PACS numbers: 03.65.Ud, 03.67.Hk 


Optical path entanglement - entanglement between 
several optical paths sharing a single photon - is one 
of the simplest forms of entanglement to produce. It 
is also a promising resource for long-distance quantum 
communication where the direct transmission of photons 
through an optical fiber is limited by loss. In this con¬ 
text, loss can be overcome by using quantum repeaters, 
which require the creation and storage of entanglement in 
small-distance links and subsequent entanglement swap¬ 
ping operations between the links. Among the different 
quantum repeater schemes, those using path-entangled 
states | 1 )a| 0 )b + | 0 )a| 1 )p 5 where a single photon is de¬ 
localized into two nodes A and B are appealing - they 
require fewer resources and are less sensitive to mem¬ 
ory and detector efficiencies compared to repeater ar¬ 
chitectures based e.g. on polarization entanglement [T]. 
Many ingredients composing these networks have been 
experimentally demonstrated, including path entangle¬ 
ment based teleportation [2] , entanglement swapping [3] , 
purification [4] , quantum storage m and an elementary 
network link [7]. 

A natural question is how this body of work could serve 
to extend known point-to-point quantum repeaters to 
richer geometries for quantum networks? FIG. [^presents 
a possible solution: A single photon incident on a multi- 
port coupler generates entanglement over N output paths 
(see FIG. [^), due to its non-classical nature [8]. The 
small network created in this way can be entangled with 
other, potentially distant, networks via entanglement 
swapping operations using 50/50 beam-splitters and sin¬ 
gle photon detectors - a single detection is then enough 
to entangle the remaining 2N-2 nodes (see FIG. Ef)- 
Such 2D networks could open up new perspectives for 
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FIG. 1: Proposal to build up 2D networks over long distances, 
a) Networks made with neighboring nodes are made with N- 
path entangled states, b) These local networks can be con¬ 
nected remotely by means of entanglement swapping opera¬ 
tions resulting in a large scale network. 


multi-user quantum information processing including se¬ 
cret sharing [9] or secure multi-party quantum computa¬ 
tion m as well as for experiments simulating quantum 
many-body system dynamics m- 

A central challenge, however, is to find an efficient, yet 
trustworthy, way to certify the functioning of these net¬ 
works, i.e. how to characterize path entanglement in a 
distributed scenario using only local measurements. One 
might do this by using several copies of path-entangled 
states, as is the case for standard quantum repeater 
schemes [12], however, doing so is resource demanding 
and addresses a restrictive class of applications - those 
accepting post-selection. State tomography has also 
been realized m to characterize two-path entangled 
states but the exponential increase in measurements 
with the number of subsystems makes the tomographic 
approach impractical for detecting the entanglement 
in large multipartite systems like quantum networks. 
Recently, an entanglement witness for bipartite path 
entangled states has been proposed and demonstrated, 
that is based on a Bell inequality combined with local 
homodyne detections m US]- However, it is not clear 
how this approach can be extended to more than two 
paths as even for three parties we know of no Bell 
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inequality that can be violated for W-like states with 
measurements lying on the equator of the Bloch sphere. 

In this letter we propose an entanglement witness 
specifically developed to reveal path entanglement in dis¬ 
tributed systems. It relies on an accurate description of 
measurement operators and assumes that each path is 
described by a single mode. However, it does not require 
post-selection, nor assumptions about the photon num¬ 
ber of the measured state, hence, it reveals entanglement 
in a trustworthy manner. Moreover, it only makes use 
of local measurements and easily scales to multipartite 
systems. 

The principle of the witness is the following: N dis¬ 
tant observers share a state p describing N optical paths. 
Assuming that each path is completely described by a 
single mode of the electromagnetic field, the aim is not 
only to say whether the overall state is entangled, but 
also to check that entanglement lies in a subspace with 
at most one photon in each mode and to check that 
pp]ni<i is genuinely entangled. The subscript i labels 

i 

the observer 1 , 2 ,....^" and rii is a non negative integer 
describing the photon number in the optical path i. To 
do this, each observer uses a measurement combining a 
small displacement operation and a single photon de¬ 
tector, a measurement initially proposed in Refs. 

[TH] and demonstrated in m- In the qubit subspace 
{|0),|1)}, the POVM elements corresponding to click 
and no-click events of such a measurement can be seen 
as non-extremal projective measurements on the Bloch 
sphere whose direction depends on the amplitude and 
phase of the displacement [20]. In other words, if one 
considers non-photon-number-resolving (NPNR) detec¬ 
tors with a quantum efficiency p and a small displace¬ 
ment D{a) = operating on the mode i, the 

corresponding observable is given by 




if one assigns the outcome +1 when the detector does 
not click and —1 when it clicks. If the measured state 
belongs to the subspace with at most one photon and 
with p = 1, 00 (the superscript is omitted when p = 1) 
corresponds to the Pauli matrix i.e. a qubit mea¬ 
surement along the z direction. Similarly, for a = 1 
and a = (uI and are a good approximation to qubit 

measurements along x and y, respectively. We use this 
analogy to build up a fidelity-based entanglement wit¬ 
ness of the form Zn = N{2^\Wn){Wn\ — 1), where 
~ •••5 1^5 •••Oat) refers to the state in¬ 

volving N modes sharing a single photon. We approxi¬ 


mate this expression by the operator 
N 

Zn = ^ (AT - 2m)(t7®™ ® 1®^-™ 

m=l 

N-2 

+2 y] (uf ” 0 0 0 (U„ + Uia 0 dTja) 

m=0 

+sym. (2) 

which only involves measurements of the form Q . Z 

t^N-m fQj. g measurement in which the first m 

paths are measured with the observable oq and the N — m 
remaining ones are traced out. ”sym.” indicates that we 
add terms corresponding to permutations over all paths. 

To make our witness suitable for experiments, we focus 
on the case where the displacements are phase averaged 
so that the relative phase of displacements is random 
but the phase of each displacement which respect to the 
state on which it operates is well controlled. Under this 
assumption, the statistics on outcomes obtained by mea¬ 
suring m paths with is the same for any cj). Hence, 

our witness reduces to 

N 

m=l 

N-2 

+4 ® ® cr„ 0 (tT„ 

m=0 

+sym) (3) 

where 0 is averaged out. In order to detect entanglement 
with Zjv, it is sufiisant to compare its value to the 
maximum value ^ d(t)TT[ZNp] that it can 

take if the projection of the state p in the {0,1} subspace 
Pp]ni<i a positive partial transposition (PPT) with 

respect to a single party. Indeed, the observation of 
a value oi Zjsi larger than implies by the Peres 

criterion |2Tl [22| that the measured state is entangled 
and that the entanglement lies in the qubit subspace. 
Since finding constitutes a linear optimization 

problem with semidefinite positive constraints, it can be 
computed efficiently (see Supplemental Material). Simi¬ 
larly, comparing the value of to the maximum 

value of z^^^g further optimized over all possible PPTs, 
reveals genuine multipartite entanglement. 

As an example, consider the value that the witness 
would take, zw, in a scenario without loss and involving a 
state Wat in which N optical paths share a single photon. 
We can compare this to the value that would be 

achieved without genuine entanglement in the {|0), |1)} 
subspace. We show in the Supplemental Material that 

= 2^+3^|a|2e-2|«l^ (4) 

which is positive for all N. The proposed witness thus 
has the capability to reveal genuine entanglement of Wat 
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FIG. 2: Three different set-ups used to test the proposed en¬ 
tanglement witness for two and three parties: a) The heralded 
state can be tuned from maximally entangled to separable 
by the half-wave-plate (HWP) Xe before the first polarizing 
beam-splitter (PBS). The local oscillator is introduced at the 
other port of the PBS such that in each arm, the coherent 
state and the single photon have orthogonal polarization. The 
displacement operation is performed by rotating the HWPs 
at Xdi and Xd2- b) The single photon and coherent state are 
input earlier in the set-up with orthogonal polarizations. The 
input loss can be varied to study the robustness of the witness. 
This set-up can be easily modified, by adding a 30/70 beam¬ 
splitter and another (dashed) arm, allowing us to herald and 
probe a tripartite W-state. 


states for any path number. In practice, the value of a 
is optimized to make the difference zvp — as large 
as possible. 


photons onto beam-splitters, see FIG. The photons 
are prepared using a heralded single photon source 
(HSPS) based on a bulk PPLN nonlinear crystal [23] . 
The crystal is pumped by a pulsed laser at 532 nm in the 
ps regime with a repetition rate of 430 MHz producing 
non-degenerate photons at 1550 nm and 810 nm via 
spontaneous parametric down conversion. The telecom 
photon is filtered down to 200 pm and subsequently 
detected by an InGaAs single photon avalanche diodes 
(SPAD), producing pure heralded single photons at 
810 nm - the purity is verified by measuring the second 
order autocorrelation function ^^(0) [24]. To ensure a 
high fidelity entangled state, the pair creation prob¬ 
ability per pulse is limited to 10“^, to minimize the 
effect of double pairs, the photons are coupled with 90 % 
efficiency m and the overall system transmission is 
optimized. We herald single-photon states at a rate of 
^ 8 kHz. 

The measurements are performed by combining 
displacement operations and single photon detection. 
The local oscillator for the displacements is generated in 
a similar PPLN nonlinear crystal pumped by the same 
532 nm pulsed laser as well as a 1550 nm telecom CW 
laser - this ensures a high degree of indistinguishability 
between the HSPS and the local oscillator, which is 
confirmed by measuring a Hong-Ou-Mandel interference 
dip between the two sources, where the visibility is 
only limited by the statistics of the two sources [24] . 
Custom gated Silicon SPADs are then used to detect 
the photons at 810 nm m- The detectors operate 
at 50% efficiency and have a dark-count probability 
of 10“^ per gate, for a gate width of approximately 2.3 ns. 

To determine the value of the witness, which reduces 
to 


When the measured state is not entirely contained 
in the {|0) , |1)} subspace, contributions from higher 
photon numbers can increase the value To get 

a valid bound in this regime, we used autocorrelation 
measurements in each mode. They give a bound on 
the probability of having more than one photon in each 
path (pc*^ denotes this bound for mode i) and avoid 
making assumptions about the photon number. The 
computation of is then slightly modified to take the 

value of pc^ into account (see Supplemental Material). 
Importantly, the autocorrelation measurements are 
performed locally with a beam-splitter and two photon 
detectors. Overall, the number of measurements required 
to reveal genuine entanglement between N paths scales 
quadratically (^ + ^ + 1)5 which shows a much more 
favorable scaling compared to the exponential scaling of 
tomographic approaches. 

We now report on a series of experiments demonstrat¬ 
ing the feasibility of our witness. We prepare entangled 
networks made with 2 or 3 paths by sending single 


^2 = (2fl7„ 0 - «To 0 00) (5) 

in the bipartite case, we measure click/no-click (c/0) 
events in the two paths and calculate the correspond¬ 
ing probabilities, Pqo, Pqc, ^cO, Pcc^ as well as the bounds 
on the probabilities for having more than one photon 
in each path, and p ^^. The correlators of the form 

{oa' in (§ then correspond to PooT Pcc Pqc PcO^ 

for a' = (a or 0. We first block the single photon from go¬ 
ing to the set-up and apply the displacement operators in 
both arms, validating that \a\ corresponds to the desired 
value. Experimentally, \a\ is such that Pc ^ Pq locally, 
see Supplemental Material. Secondly we allow the sin¬ 
gle photon to go to the set-up and record the correlators 
with, and without, the displacement operations. An au¬ 
tomated series of measurements is performed, integrating 
over 1 s for each setting, and is repeated as many times 
as needed to have good statistics. The values for p^^^ and 
Pc are dominated by detector noise due to operating the 
detectors at such high efficiencies so as to maximize the 
global efficiency of the measurements. These values are 
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FIG. 3: Observed value for the bipartite entanglement witness 
(relatively to the PPT bound) as a function of the beamsplit¬ 
ter ratio. Concretely, the half waveplate \e in Fig 2a) is 
rotated which changes the state from a maximally entangled 
to a separable state (50/50 - 0/100 splitting ratio, respec¬ 
tively). The blue band is obtained from a theoretical model, 
taking into account the set-up’s global transmission, the char¬ 
acteristics of sources and gated detectors and the value for 
\a\ ~ 0.83 in the displacement operations. 


used to determine the observed value of labelled 
and the maximum value that would be obtained if 
the projection of the measured state in the {|0), |1)} sub¬ 
space has a positive partial transpose (see Supplemental 
Material). 

To test the bipartite witness as a function of the 
amount of entanglement, the single photon and local 
oscillator are combined at different ports of a polarizing 
beam-splitter (PBS) ensuring that they leave in the 
same spatial mode with orthogonal polarizations, see 
FIG. 1^. A half-wave-plate (HWP) Xe placed in the 
single photon input arm is used to adjust the splitting 
ratio in the two output modes and the subsequent am¬ 
plitudes for the entangled state, are performed via a 
rotation of the wave-plates Xdi^ Xd 2 (< 1 degree) before 
the final PBSs. The amplitude of the displacement 
\a\ 0.83 is set to maximize z^^^ — z^^^. FIG. 

shows the result as a function of the beam-splitter ratio, 
from maximally entangled (50/50) to a separable state 
(0/100). The shaded line is obtained from a theoretical 
model with independently measured system parameters, 
with the associated uncertainty (see Supplemental 
Material). The theory and experimental results are in 
excellent agreement and prove that the proposed wit¬ 
ness can reveal even very small amounts of entanglement. 

To prove the robustness of this witness against loss, 
and demonstrate the scalability, we introduce a different 
experimental configuration, FIG. |^, with only a 50/50 
beam-splitter, to generate maximally entangled states, 
and the single photon and local oscillator are combined 
earlier in the set-up. We can then introduce loss to the 
input state, thus increasing the mixedness of the state. 
FIG. shows the value of our witness of entanglement as 



Total Transmission 

FIG. 4: Observed value for the bipartite entanglement witness 
(relatively to the PPT bound) as a function of loss. The to¬ 
tal transmission consists of the photon coupling, transmission 
through the system, and detector efficiencies. 


we increase loss. The starting point has a slightly larger 
value than in FIG. due to a slightly better transmis¬ 
sion, n.b. the maximum transmission of >30 % includes 
photon coupling, transmission and detection efficiency. 
Here we see that the witness is capable of revealing en¬ 
tanglement even in the case of high loss, or similarly, for 
low detection efficiency. 

Finally, by adding a 30/70 beam-splitter and another 
arm, dashed line in FIG. It, we herald tripartite 
states. If we assume perfect transmission and detectors 
with unit efficiency, we expect a maximum value for 
— z^^^ ^ 7.63, where a value greater than zero 
indicates the presence of genuine entanglement. By 
applying our model, similarly to the bipartite case, again 
with \a\ ^ 0.83, but with a total transmissions in each 
arm of 0.19 ± 0.002, we expect to find a theoretical value 
of 0.99 (see Supplemental Material). We found a value 
of 0.99T0.10 that agrees with our model and shows a 
clear violation, thus revealing genuine tripartite path 
entanglement. 

In conclusion, we have shown an entanglement witness 
suited for path entangled states that is robust and 
scalable, providing the means for the characterization 
of genuine multipartite entanglement distributed over 
complex quantum networks. The co-propagation of the 
local oscillator with the path entangled state overcomes 
the potential problem of distributing a phase reference, 
which also has the added advantage that it could 
be exploited for stabilisation and synchronisation of 
distributed networks. Interestingly, our witness provides 
a trustworthy means to reveal entanglement, without 
the need to make assumption about the number of 
photons in each path. A possible extension would be to 
make it fully device independent through the violation 
of a Bell inequality, which would require higher overall 
efficiencies |23|27]. 
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Appendix 



A. PPT bound for an arbitrary number of qubits 

In this first paragraph, we give details on the result 
presented in Eq. (4) of the main text. Our witness has 
the following form 

N 

m=l 

N-2 

+4 ® ® tJa ® 

m=0 

( 6 ) 

The question we address is what is the maximal value 
Tr[^ dcjyZNp] that Zn can take if the projection of 
the state p in the {0,1} subspace Pf|ni<i stays positive 

i 

under partial transposition with respect to some bipar¬ 
tition. In this section, we consider the case where each 
mode is described by qubits p = p{^ni<ii i-C- they con- 

i 

tain at most one photon each. This condition is relaxed 
in the next section. The threshold value that we look for 
can be obtained from the following optimization 

1 

= inax Tr[— / #^ArPnni<i] s.t. 
Pnni<i Jo i ~ 

(^) Pn 

i 

(ii) Tr/)Q„,<i = 1, 

i 

(Hi) Pfjn-^i ^ 0- 

The two first conditions insure that pf^rii<i is a physical 

state. The last one demands that the state remains 
positive under partial transposition with respect to the 
first path. This is a linear optimization with semidefinite 
positive constraints which can be efficiently calculated 
numerically. If a physical state made with qubits leads 
to a larger value than we conclude that the 

condition (iii) does not hold, i.e. its partial transposition 
over the first path has at least one negative eigenvalue 


FIG. 5: Value that the witness takes with a Ww state zw rela¬ 
tively to the maximal value with states that are not genuinely 
entangled 2 :^^^ as a function of the number of optical paths. 
The absolute value of the difference — does not matter 
a,s Zn could be multiplied by any values. However, a positive 
difference — 2 ;^^^ > 0 reveals genuine entanglement. 


and is thus entangled. If we further optimize over all 
possible bipartitions, we obtain a bound from 

which one can witness genuine entanglement between 
the A optical paths. 


Consider for concreteness the case where the tested 
state is the Wn state. We can easily show that zw = 
C \Wn) (Wiv|] is given by 

zw = {2^-l)N+ 2^+^{N (7) 

X (|a|2(4e-l“l'-l)-l). 

Furthermore, a basic measurement of the photon number 
distribution in each path would show that the probabil¬ 
ity of having zero photon in all paths or more than one 
photon per path is null (see below for the realization of 
such measurement). Together with the condition (iii), 
this leads to conditions on the coherence terms of the 
tested state, e.g. 

|(10...0|p|01...0)|2<poo...oPii...o = 0. (8) 


The conditions (i), (ii) (iii) together with the previously 
mentioned constraints allows us to get an analytical ex¬ 
pression for the PPT bound with respect to the first path 


(9) 

More generally, if we consider the entanglement between 
the first m modes and the other (N—m) modes by chang¬ 
ing the condition over the partial transposition, we obtain 




( 10 ) 


which gives the threshold for a m vs. {N — m) bisepara- 
ble state. This bound is identical for all such bipartitions 
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because our witness and the conditions (i)-(ii) are invari¬ 
ant under exchange of parties. The threshold for genuine 
multipartite entanglement can be deduced by minimizing 
zw — over m. This leads to 

= (11) 

which is always positive for non-zero \a\. FIG. 1 is the 
result of the optimization of the previous expression over 
a as a function of the number of parties. This shows 
that our witness can detect genuine entanglement of IFw 
states for any number of parties. 


B. PPT bound for qudits 

We consider the general case where no assumption is 
made about the number of photons in each path (no 
qubits). If the measured state is described by a single 
mode, it can be written as 


PC]ni<l Pcoh 

P = \ ' t 

Pcoh P[J'ni>2 


( 12 ) 


where is the probability to have strictly more 

than one photon in the path 1. Since p (y Hi >2^ < 
EZi we have 


P (Ijn* > 2j < 2(p(E +pZ^ + ... +pi^^). (15) 

The previous formula shows how to bound the probability 
for being outside the qubit subspace {|0), with 

local measurements involving a 50/50 beam-splitter and 
two photon detectors. From Eqs. and ( p!^ , we can 
deduce the PPT value in the general case of qudits from 
the following optimization 


+2^,lg(p(l) +p(2) + ...+pm) S.t. (16) 
(0 Pnni<2 > 0, 

i 

(ii) Tr/)Q „.<2 < 1, 

i 


where Pf^rii<i denotes the block with at most one photon 
per partie as before, p\jni >2 is the block in which at least 

i 

one party has more than one photon and pcoh denotes the 
coherence terms between these two blocks. Importantly, 
the only term in pcoh leading to a non-zero contribution 
to our witness are those between pf^rii<i p\jni= 2 - 

Given the maximal algebraic value z^ig that can take, 
we have 



/ 

i. 

p27T 

/ 

Pr\ 

i 

Pcoh 0 \ 

\ 

zn < Tr 

/ 

Jo 


J 

Pcoh 

P\Jni=2 0 



\ 

\ 

0 

0 

0 

/ 


T ZoigP ^L_J ~ 

(13) 

where p ^|J > 2^ is the probability to have more than 

one photon in at least one path. This probability can 
be upper bounded using a 50/50 beam-splitter followed 
by two detectors in each path. Indeed, let us consider 
the reduced state describing, say, the first path pi. Its 
diagonal elements in the Fock basis can be written as 
Sn>o 1^) (^1 • The probability pi^^ that this state leads 
to a two-fold coincidence after the beam-splitter is given 

by 

i-i" = E |r( 2 ”- - i)rf> i (») 

n>2 n>2 


Note first that the optimization is performed over the set 
of physical states (condition (i)) having two photons at 
maximum in each path because the only coherence terms 
having a non-zero contribution lie in this subspace. To 
keep the optimization general, the trace of the state is 
required to be smaller than, or equal to, one (condition 
(ii)). The last condition ensures that the projection of 
the state with at most one photon in each path has a 
positive partial transpose when the transposition is taken 
over the first path. By changing the PPT condition (iii) 
to > 0 where Tm denotes partial transposition 

with respect to the bipartition m, we get a PPT bound 
for this bipartition The value of the PPT bound 

we are interested in is = max^ Importantly, 

the difference between and is made larger 

by using the results of the joint probability to have 
clicks/no-clicks without displacement and adding the 
corresponding constraints in the optimization procedure, 
cf below. Overall, the number of required measurements 
(1 + C 2 ) + = A^^/2 + N/2 + 1 has a quadratic scaling 

in the number of paths and is thus suited to prove 
genuine entanglement in multiple-path states. 


C. Witnessing path entanglement with non-unit 
efficiency detectors 

Note that so far we have assumed that the photon de¬ 
tectors have unit efficiencies. How can we use this witness 
in practice, when non-unit efficiency detectors are used? 
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First, note that a detector with efficiency r] can be seen 
as a unit-efficiency detector preceded by a beam-splitter 
with a transmission efficiency r]. Let t/bs be the unitary 
corresponding to this beam-splitter. Let also D{a) be the 
unitary associated to the displacement operation with ar¬ 
gument a. We have 

D{a^)U^, = UbsD{a) (17) 

meaning that the inefficiency of the detector can be 
modeled as a beam-splitter operating before the dis¬ 
placement operation provided that the amplitude of 
the displacement is reduced by Let us consider 

the configuration where the loss operates before the 
displacement. The optimizations performed so far allows 
one to conclude that the state after the loss is entangled. 
This implies that the state was already entangled before 
the loss as a separable state remains separable under 
loss. Our witness thus proves entanglement when 
non-unit efficiency detectors are used provided that 
the displacement accounts for the reduced detection 
efficiency. 


D. Bipartite case 

We now give an explicit expression of the PPT bound 
for the bipartite case. For = 2, we take the following 
witness 

Z 2 = (2(t7„ <S)€a-€o<S) 0o) (18) 

Note that we have removed a factor 2 from the definition 
To get the PPT bound first note that 

tne maximal (algebraic) value for Z 2 is 3 and since 
P {ni > 2 U ni > 2) < 2{p^c^ we have 


^max 

^ppt 


< max Tr ( T / d<t)Z2pm<2nn2<2 

Pni< 2 nn 2<2 \ Z7r Jq 




s.t. 


f- Pni<2P|n2<2 ^ O 5 
2 . Tr(p^^< 2 p|n 2 < 2 ) < I 5 

^ni<inn2<l — 

4. 1 -p{ni > 2[jn2 > 2) = Tr(p^^<i 

When oo ^ is measured, one accesses the joint prob¬ 
abilities of having click/no-click without displacement. 
Let Pij be the joint probability to have outcomes i and 
j respectively, i,jf = 0 for no-click and = c for 
click. The Pijs provide upper bounds of the diagonal 
terms of the measured state poo = Poi ^ Poc^ 

Pio < Pco and pu < Pec- Similarly, from the local 


probabilities of having more than one photon, we have 
P02 < ‘2p^c\ P20 < Pl2 < P 21 < 

and P 22 < 2p^\ Together with the condition 1. and 
3. (implying e.g. |(10|p„^<2nn2<2|01)|^ < PoiPio and 
|(10|/9„i<2nn2<2|01)P <PooPii respectively), these con- 
straints provide the following upper bound on in 

the regime a > 0.45 

< (^2 (-l + 2e-W^)'-l)Poo 

+ (2 (-1 + 2e-l“l') (-1 + 2e-l“l' jap) + l) (Poc + Pco) 

+ 2(2 (-1 -t - 1^ 

+ 2 (2 (-1 + 2e-l“l") (-1 + |a|^) + d) { p ^^^ + p^^)) 

-I- (^-\/pP^\a\‘^ + \/PqoPcc 

+ (^\fjppcc + \fjppcc + \fPpPj |ap)- 

(19) 

This analytical bound is used and compare to the ob¬ 
served value of our witness. If the latter is larger than 
the former, we can safely conclude that the tested state 
is entangled. 


E. Theoretical model 

The aim of this section is to detail the model that has 
been used to reproduce the experimental results (blue 
lines in FIG. 3 and 4 of the main text). The starting 
point is the pair production in two different modes by 
means of a spontaneous down conversion source and sub¬ 
sequent photon detection of one of two modes. The state 
conditioned to a click on the heralding detector is given 
by [ 20 ] 


Ph = 


1 - RjT^ r 


THI-Rl) 


Pt\i\n = 



1 - T2 (_ _ 


l-RlTp\ 


( 20 ) 


i.e. a difference between two thermal states pth{P) 
with n mean photons. Tg = tanh(p), g being the 
squeezing parameter, and a stand for bosonic op¬ 

erators and Rh = \/l — Ph where ijh is the heralding 
efficiency. Importantly, a thermal state can be writ¬ 
ten as a mixture of coherent states jq) . Concretely, 

Pth{n) = f i'll with P”( 7 ) = Ve-4J-, 

The correlators that we want to calculate can thus be 
obtained from the behavior of a coherent state. A beam¬ 
splitter splits a coherent state into two coherent states, 
i.e. I 7 ) ^ |V^ 7 )q IvTq)^, where T and R are, respec¬ 
tively, the transmittivity and reflectivity. A displacement 














D{a) on a coherent state I 7 ) gives another coherent state 
with mean photon number a\‘^. Together with 

1 ^) = 1 ^x 7 ), ( 21 ) 


we easily obtain the probability to get no click in both 
sides from a thermal state pth{n) knowing the amplitude 
of the local displacement a 


-2\a\^ + j^^{VR+VTf\a\^ 


TDOi _ 

^00 — 


1 + nr]tr] 


( 22 ) 


Finally, the value taken by our entanglement witness 
is obtained from 2Tr(p/i©-ci 0 Oci) ~ ^ ®’o)- 

An independent characterization of the source to¬ 
gether with the measurements of efficiencies ph and 
PtP shows a very good agreement between our model 
and the results of the experiments (see FIG. 3 and 4 
of the main text). The error associated to both the 
theoretical models in FIG. 3 and 4 of the main text 
take into account the error on the value of a plus the 
error in the transmission of the set-up for FIG. 3 and in 
the beam-splitter ratio that prepares the state for FIG. 4. 


where ptV is the overall efficiency (from the source to the 
detector, including the detection efficiency). Following 
the same line of thought for P^, Pq^, P^, we find 


Tr(pth(^)^^a G (SJa) =1 + 4 


1 + nptP 


p l-\-r]ir]nR p l-\-r]ir]nT 

-z-2- 

1 + 1 + PtP'riT 


F. Tripartite case 


For = 3, the expression of the witness Z 3 is 


From this last expression, we deduce the correlator for 
the state ( 20 ) 


1 - RlT^ 

Tr {ph€a ® Ua) = ^ X 


THi-Rl) 


Tr Pth 


rp2 

9 


1 - T2 

9 , 


OTq/ yy Oq; 


1 - 7 ( 

® Tr Ah 




' \ 
l-RlT^ 


(23) 


J a . ^ 


Zs = 3oo G ^Jq) 

+4 (1 + Oo) <Ua ^ <^a + sym^ ^^ 24 ) 


where sym stands for the terms obtained by all possible 
permutations of modes 1, 2, and 3. Following the same 
line of thought as for the bipartite case, we can show that 
for any a > 0.67 


^ppt 


)) X (Pooc 


< (^—3 + 24 (“1 + 2 e ^ ^ -Pooo + (s + 16 (~1 + 2 e ^ (~^ ^ {Rooc + -foco + Rcoo) 

+ 4(l + 4(-l + e-l“l+|^) (-3 + 4e-l“l'x (pi^^+p7+p7) 


+ 64|a|2e-2|“l' (jap( + VP+F+ + ■\/p+F+ + VP+F+) 

+ \a\^ (l + lap) (^^4+++ 

+ max(-\/Poco-Fboc + \/PoooPcOc + VPoooPcco, \J -fcoo-Fboc + \J -Pooo-Fbcc + \lPoooPcco, 
\/a)c0-Pc00 4“ RooqRcOc 4“ \/-Tooo Rqcc ) ) 

+ 33(2 +pR +p 7 ) - max (0,2 {^R - 1,2 {pR - 1,2 [pR “ ^) ) 


Observation of a value larger than allows one to G. Experiment 

conclude that the measured state is genuinely entangled. 


FIG. 2 shows the experimental scheme for the 
generation of heralded single photons and the local 
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FIG. 6: The experimental set-up for the heralded single pho¬ 
ton source and the local oscillator. See text for details. 


oscillator for the displacement operations. Most of the 
experimental parameters are explained in the main text. 
We have two independent sources based on PPLN non¬ 
linear crystals, with type-0 quasi-phasematching. The 
heralded single photon source (HSPS) uses a tunable 
filter (F) with a 200 pm bandwidth that has around 
50% transmission for the central (telecom) wavelength. 
This is sufficiently narrow so as to herald the 810 nm 
photons in pure states M- By stimulating the DFG 
process in the second non-linear crystal with a coherent 
CW laser at the same wavelength of the HSPS’s filter, 
energy conservation dictates that both the heralded 
single photon and the local oscillator (LO) state are 
indistinguishable. A delay arm (not shown) allows for 
the sychronization of the two sources. 


To determine the value of \a\ which maximizes 
the probabilities Poo, ^Oc, ^co, ^cc, as well 

as and are measured. By guessing that they are 
obtained from a W state, these probabilities allow one 
to estimate the overall efficiencies and hence the value 
that the witness would take for any a. Maximizing the 
difference with (19) results in \a\ = 0.83. We emphasize 
that in practice, with non-unit efficiency detectors, this 
corresponds to a displacement with amplitude 0.83/y^. 
This means that the detection efficiency does not need 
to be known. The values of |(a| are locally set to produce 
click/no-click with the same probability independently 
of the detection efficiency. 


The stability of the mean photon number for the LO 
states is realized by a feedback loop, where a large frac¬ 
tion of the generated LO is sent, via a polarizing beam¬ 
splitter (PBS) to a power meter (PM) that then acts on 
a Pockels cell. The typical integration times, per point in 
FIG. 3 & 4 in the main text, are around 1.5 hours. The 
error bars are not Poissonian as they are dominated by 
fluctuations in the system, in particular, the generation 
of the LO, which depend on the feedback loop. 

If we consider the measurement for the maximally en¬ 
tangled bipartite case as an example, the probability 
to have more than one photon in each arm are mea- 
1.0 


0.8 

-M 

^ 0.6 

03 

_Q 

O 0.4 

L- 

Q_ 

0.2 

0.0 


FIG. 7: Measured probabilities of having click/no-click in the 
two arms of the bipartite experiment (FIG. 4 of the main 
text with the maximal transmission efficiency (> 0.3)). Red 
columns give the probabilities measured without the displace¬ 
ment operation and the blue columns correspond to the prob¬ 
abilities measured when the displacement operator was ap¬ 
plied (|(a| = 0.83 ± 0.01). 



sured to be p^^^ = 10“^ and p^^^ = 10“^. After mea¬ 
suring the different click/no-click events we can deter¬ 
mine their probabilities - see the results presented in 
FIG. 0 These lead to = -0.002, = -0.315, 

i.e. i^P — = 0.313 > 0. As the observed value 

was larger than ^^p^p^, we can conclude that the measured 
state is genuinely entangled. 
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